Two-Step Restoration of SU(2) Symmetry in a Frustrated Ring Exchange Magnet 
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We demonstrate the existence of a spin-nematic, moment-free phase in a quantum four-spin ring 
exchange model on the square lattice. This unusual quantum state is created by the interplay of 
frustration and quantum fluctuations which lead to a partial restoration of SU (2) symmetry when 
going from a four-sublattice orthogonal biaxial Neel order to this exotic uniaxial magnet. A further 
increase of frustration drives a transition to a fully gapped SU (2) symmetric valence bond crystal. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.40.Cx 



Broken symmetries are one of the central paradigms of 
magnetic ordering, and most antiferromagnetic systems 
are in Neel phases at low temperatures, characterized by 
a vectorial order parameter: their sublattice magnetic 
moment. This order parameter breaks the rotational 
SU{2) symmetry of the Haniiltonian and is accompanied 
by gapless excitations, the Goldstone modes of the broken 
symmetry Q . In low-dimensional systems of Heisenberg 
spins frustrating couplings can drive transitions to gap- 
ful SU{2) symmetric quantum states, where the building 
blocks are local singlets (in the simplest examples, pairs 
of spin-1/2 in short-range singlet states). These quan- 
tum gapped phases may have long-ranged singlet order 
and break spatial symmetries of the lattice, called va- 
lence bond crystals (VBC) in the following. An even 
more exotic groundstate, a coherent superposition of all 
lattice-coverings by local singlets, a state known as a res- 
onating valence bond (RVB) spin liquid, could also be 
found 2]. 

Quantum Phase Transitions from Neel ordered phases 
to quantum gapped phases have been studied for a long 
time as prototypical examples of quantum phase transi- 
tions. The nature of the symmetry breaking Neel phase 
plays an important role in these scenarios and seems de- 
terminant as to whether the adjacent gapped phase will 
be a VBC or an RVB phase jSj. It was recently shown 
that the transition from the standard coUinear (tt, tt) Neel 
state to a Valence Bond Crystal phase can actually be 
an exotic quantum critical point with deconfined excita- 
tions 3. 

The well-known (tt, tt) Neel state is a uniaxial magnet, 
with two gapless Goldstone modes, in which SU{2) is 
partially broken to U{1). More complete SU{2) break- 
ing schemes do exist, for example in noncoUinear magnets 
with more than two ferromagnetic sublattices or more 
generally in helicoidal antiferromagnets. In these sys- 
tems, the order parameter can be described as biaxial 
(or as a top), the SU{2) symmetry is completely broken, 
and there are three Goldstone modes. Chandra and Cole- 
man suggested that in such situations the restoration of 



the full SU{2) symmetry due to the interplay of quantum 
fluctuations and frustration could possibly occur in two 
steps: from a biaxial magnet via an intermediate uniax- 
ial spin nematic magnet - still with gapless excitations 

- to a fully gapped paramagnetic phase without SU (2) 
symmetry breaking j^, |^ ■ This speculated spin-nematic 
phase - first introduced by Andreev and Grishchuk j3| 

- has no net magnetic moment, but nevertheless breaks 
the SU{2) symmetry, as the individual spins - albeit dis- 
ordered - remain correlated in a plane. Up to now this 
conjectured two-step scenario misses a concrete realiza- 
tion, as it has not been found in the models originally 
proposed, e.g. the Ji — J3 model on the square lattice or 
various models on the kagome lattice. 

In this Letter we present the phase diagram of a frus- 
trated four-spin ring exchange model for S ~ 1/2 on 
the square lattice. We show that this model has a four- 
sublattice orthogonal Neel groundstate with a biaxial or- 
der parameter. This groundstate maximizes the square 
of the vectorial chirality, and appears as a natural "un- 
frustrated" starting point in the following. Increasing 
the frustration by a change in the nearest neighbor an- 
tiferromagnetic coupling, provides a first realization of 
the mechanism speculated by Chandra and Coleman |^ . 
Based on the analysis of spin-resolved spectra, the transi- 
tion from a biaxial to a uniaxial magnet - accompanied by 
the reduction from three Goldstone modes in the biaxial 
magnet to only two Goldstone modes in the uniaxial one 

- is highlighted. Upon a further increase in the frustra- 
tion all spin excitations become gapped and the system 
enters a (staggered) Valence Bond Crystal phase. 

The model Hamiltonian with spin 5 = 1/2 reads: 

H^J J2 S,-S,+K J2 {P^...i + P^'.i) ■ (1) 

<i,j> [i,j,k,l] 

Pi,,,i is the cyclic permutation of the 4 spins sitting on 
a square plaquette, and the two sums run respectively 
over the nearest neighbor bonds and the plaquettes of 
the square lattice. We parametrize the couplings by J = 
cos 9 and K = s\\\9. The explicit expression of the 4-spin 



2 



K 

I Spin 
orthogonal Nematic 




I 



FIG. 1: (Color online) Schematic phase diagram of the cyclic 
four-spin exchange Hamihonian Q on the square lattice as a 
function of the parameter 6. 



permutation operators in terms of spin-1/2 operators can 
be found in Ref. 

Interest in the present model has arisen 15 years ago 
in an early attempt to describe the magnetism of the 
cuprates[9|, [lOj. It has been shown very recently to be 
relevant to describe the high-energy spin waves disper- 
sion of La2Cu04 In this material, the four-spin ring 
exchange term is significant, but still rather small com- 
pared to the usual Heisenberg term jl^j. Here we study 
this model in a broader range of parameters. Recent 
studies of the Hamiltonian |Q on a two-leg ladder have 
provided evidence for a number of unconventional phases 
|l3l . Of particular relevance for the present work is 
the appearance of a short range ordered, gapped phase 
with dominant vector chirality correlations in a large re- 
gion of the phase diagram. 

Our results on the overall phase diagram obtained 
by large scale exact diagonalizations of systems up to 
TV = 40 spins are summarized in Fig. ^ Adjacent to the 
well established Neel phase around J — 1. i^T = 0, we 
find two VBC phases, a staggered dimer phase for posi- 
tive K and a columnar dimer phase for negative K. For 
negative K we further find a large ferromagnetic region. 
These phases ioi K < are analogous to those observed 
in a related XF-like model, amenable to Quantum Monte 
Carlo simulations ISj]. The emphasis of the present pa- 
per is on three of the phases found for K > 0: the copla- 
nar orthogonal antiferromagnet, the spin- nematic state 
and the staggered dimer VBC. A detailed discussion of 
the whole phase diagram will be presented in a forthcom- 
ing article 16]. 

The orthogonal four-sublattice AFM — Let us start 
the discussion with a state which has a simple classi- 
cal interpretation. The cyclic exchange term K contains 
both four-spin couplings and two-spin exchange terms. 
For J = — 2 and if = 1 (0 w O.SStt), the effective near- 
est neighbor exchange coupling is vanishing . Classical 
and semiclassical reasonings then predict an orthogonal 
four-sublattice AFM [l|l3|, as shown in Fig.EJi). It could 
be remarked that this state minimizes the SU (2) and lat- 
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TABLE I: Irreducible representations of the square lattice 
space symmetry group appearing in the tower of states of the 
orthogonal Neel state (A) and the spin-nematic state (B) 
respectively, k denotes the momentum, 71^/2 the 90° rotation 
around a site, and a the site-based reflection along the x or y 
axis. N is the number of sites of the sample. 
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FIG. 2: (Color onhne) (a) The orthogonal Neel state. It has 
a coplanar, four-sublattice structure. The vector chirality on 
a plaquette possesses a (tt, tt) structure, (b) Tower of states in 
the orthogonal four-sublattice Neel state for J = —2, K — 1 
{9 « O.SStt). The eigenstates between the dashed lines form 
the groundstate manifold of the orthogonal Neel state. Taking 
into account exact degeneracies of some of the states we find 
25'+ 1 low- lying levels in each spin-S' sector, consistently with 
a complete SU (2)-symmetry breaking order parameter. Inset: 
finite size scaling of the spin gap, indicating a vanishing spin 
gap in the thermodynamic hmit. 

tice symmetric Hamiltonian: Ti' — — ^ ^ C'^ , where C is 
the vectorial chirality on an elementary plaquette: 

C = Si AS2 + S2AS3 + S3AS4 + S4AS1. (2) 

Interestingly this Hamiltonian is rather close to our 
model Eq. Q in the region J = — 2, K — 1. 

We now turn to the 5=1/2 fully quantum case, where 
the presence of the orthogonal Neel state has not yet been 
shown explicitly. We investigate the nature of the collec- 
tive ground-state with a spectroscopic method which re- 
lies on the analysis of symmetries and scaling of the low 
lying levels of the exact spectra of the Hamiltonian 
(called Quasi Degenerate Joint States (QDJS) in the fol- 
lowing) [13, 113 . These QDJS embody the dynamics and 
symmetry of the order parameter of the orthogonal Neel 
state. As it is a biaxial magnet, its order parameter is a 
quantum top. The leading term of the free dynamics of 
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a quantum top is proportional to the square of the to- 
tal spin of the sample S: its effective spectrum involves 
(25'+ 1) distinct eigenstates in each S sector, with eigen- 
values scaling as 5(5 + 1)/A^. 

Fig-Eb) indeed displays such a tower of low lying levels 
well separated from the other excitations. The symme- 
tries of the QDJS (displayed in Tab.QJ are those predicted 
by the ab initio symmetry analysis; three soft modes at 
(0, tt), (tTjO), (tt, tt) signal the full symmetry breaking of 
SU{2). The finite size scaling of the QDJS is regular and 
as expected, the tower of states collapse to the ground- 
state as (Inset in Fig and [ig). Long wave- 
length quantum fluctuations, estimated in a spin wave 
approach, lead to a reduction of ^ 30% of the sublattice 
magnetization the thermodynamic limit. The real-space 
spin correlations as well as the vector chirality correla- 
tions are in perfect agreement with these results. Based 
on the analysis of the exact spectra and finite size scaling 
of the orderparameters we believe that the four-sublattice 
Neel phase is stable for 0.4 tt < < 0.9 tt. 

The spin-nematic phase — Frustrating the four- 
sublattice orthogonal state by increasing J induces a 
drastic modification of the low lying spectrum of Eq. ^ 
which evolves towards the typical behavior of Fig. 
The 1/N finite size scaling of this tower of states 
proves that this phase breaks SU (2) symmetry [Inset of 
Fig. Eb)] • But the QDJS which display only one level in 
each 5 sector, embed the dynamics of a rigid rotator: the 
magnet is a uniaxial magnet, i.e. SU(2) is only broken 
down to U{1). One observes an enlargement of the spa- 
tial symmetry of the order parameter (see column (B) of 
Table nj, incompatible with a standard (tTjTt) antiferro- 
magnet, but consistent with a staggered long range order 
in the vectorial chirality |2Jl. This is confirmed by the be- 
havior of the correlations in the bond chirality (defined 
as V{i,j) = {Si A Sj)) shown in Fig. Eb)- On the other 
hand the finite size scaling of the spin-spin correlations 
points to a wiping out of the sub-lattice magnetization 
by long wave-length quantum fiuctuations. Such a state 
is therefore a p-spin-nematic state H, IS 0; character- 
ized by the absence of any sublattice magnetic moment 
(Si) = 0, and by the presence of a pseudo- vectorial order 
parameter V{i,j) ^ 0. 

The partial restoration of the SU (2) symmetry when 
going from the four-sublattice orthogonal state to the 
nematic state can be tracked by plotting the relative mo- 
tion of the different symmetry-breaking levels within the 
tower of QDJS while lowering d. The energy differences 
displayed in Fig. 0] show how all but one level for each 
spin sector evaporate once Qj-K < 0.5. Since the sym- 
metry group of the orthogonal four-sublattice antiferro- 
magnet is contained in the symmetry group of the spin- 
nematic state we might expect the transition between the 
two states to be a continous quantum phase transition, 
although this remains an open problem. 
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FIG. 3: (Color online) (a) Real space vector chirality cor- 
relations ([So A 81]"= [Si A Sj]^) for a = 40 sample in the 
spin-nematic phase sX Q = O.Stt. The black bond denotes the 
oriented reference bond. The width of the lines is propor- 
tional to the correlation strength, (b) Tower of states in the 
spin-nematic state. Inset: finite size scaling of the spin gap, 
indicating a vanishing spin gap in the thermodynamic limit. 
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FIG. 4: (Color online) The evolution of the finite size spectral 
gaps within the QDJS of the orthogonal Neel state on a = 
32 sample. The bold lines denote levels which remain in the 
QDJS of the spin-nematic state. The other levels detach from 
the QDJS as 6) < 7r/2. 



The finite size scaling of the order parameter indicates 
that the phase should at least exist in the range of pa- 
rameters 0.25 < ^/tt < 0.4. The accuracy in the determi- 
nation of the boundaries cannot be made better on the 
basis of exact diagonalizations. 

The staggered dimer VBC phase — Once the nematic 
state has been destabilized by even stronger frustration 
we find evidence for a VBC state with a staggered dimer 
structure. We consistently see an increase of the stag- 
gered dimer structure factor for all system sizes consid- 
ered. The real-space dimer correlations for an A'^ = 36 
sample are shown in Fig. [SJl)- These correlations show 
a clear staggered pattern and they converge to a finite 
value at the largest distances. Another strong argument 
in favor of a staggered dimer phase is the presence of 4 
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FIG. 5: (Color online) Staggered Valence Bond Crystal 
phase: (a) Dimer correlations in the groundstate of a 36 site 
sample at S « O.lTvr. Full, red (dashed, blue) lines denote neg- 
ative (positive) correlations. The line width is proportional 
to the correlation, (b) Low energy spectrum of a 36 sites sam- 
ple in the same phase. The required four singlets (the singlet 
at (tt, tt) marked with a star is two-fold degenerate) with the 
correct quantum numbers corresponding to a staggered VBC 
are found below the first triplet. 



singlets which will form the fourfold degenerate ground- 
state manifold in the thermodynamic limit, as shown in 
Fig. ISjs). We find two distinct singlets at momentum 
(0, 0) and a doubly degenerate singlet at (tt, tt). The sym- 
metry properties of these singlets precisely correspond to 
those expected for a staggered dimer phase. We note 
that these findings are in line with earlier work which 



proposed a staggered dimer phase in a cyclic exchange 
model on the square lattice [l^ and recent studies of a 
two leg ladder U^^. 

Conclusion — A thorough examination of the low lying 
spectrum of the ring exchange Hamiltonian on the square 
lattice revealed how a biaxial magnet may be driven to 
an uniaxial spin-nematic phase through the interplay of 
frustration and quantum fluctuations. These fluctuations 
drive the disappearance of the net magnetic moment of 
the planar orthogonal four-sublattice state, the spins dis- 
order in the spin plane, the associated Goldstone mode 
acquires a gap, but nevertheless in a a finite range of 
parameters the plane of spins remains locked. To our 
knowledge, this is the first clear demonstration of the ex- 
istence of such a phase in a two-dimensional quantum 
magnet. At last, increasing again the frustration, the 
SU{2) symmetry is completely restored, a spin gap opens 
simultaneously with the collapse in the 5 = sector of 
the four levels leading to the spatial staggered VBC. This 
is one realization of the two-step restoration of symmetry 
speculated some years ago by Chandra and Coleman p. 

In some aspects, such as spin susceptibility and ther- 
modynamic properties, the spin nematic phase is not dif- 
ferent from a standard Neel phase 01 ■ But the absence of 
long range order in ordinary spin-spin correlation func- 
tions implies that such systems - although ordered - do 



not display a Bragg peak with unpolarized neutrons. Due 
to the pseudo-vectorial nature of the order parameter, 
the NMR pattern of the nematic state is different from 
that of an ordinary Neel phase (different selection rules), 
providing a clear experimental signature of this phase. It 
might, however, be obscured by a phenomenon not ad- 
dressed in this paper, which is the sensitivity to disorder 
and impurities: disorder might be able to locally pin the 
spins in the transverse plane leading to a novel type of 
spin glass. 
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